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Abstract. This letter presents a numerical study of the evolution of an emission line cloud of initial density 10 cm~^, temper- 
ature 10"* K, and size 200 pc, being overtaken by a strong shock wave. Whereas previous simple models proposed that such a 
cloud would either be completely destroyed, or simply shrink in size, our results show a different and more complex behaviour: 
due to rapid cooling, the cloud breaks up into many small and dense fragments, which can survive for a long time. We show 
that such rapid cooling behaviour is in fact expected for a wide range of cloud and shock properties. 

This process applies to the evolution of emission line clouds being overtaken by the cocoon of a radio jet. The resulting small 
clouds would be Jeans unstable, and form stars. Our results thus give theoretical credibility to the process of jet induced star 
formation, one of the explanations for the alignment of the optical/UV and radio axis observed in high redshift radio galaxies. 

Ksy words. Galaxies: jets - Galaxies: high redshift - Galaxies: active - Galaxies: evolution - Cosmology: early Universe 
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1. Introduction 

Regions of luminous optical line and continuum emission near 
high redshift radio galaxies {z > 0.6) are often found to be 



extended along the direction of the radio axis (Chambers et al 



1987; McCarthy et al. 1987). One obvious explanationfor these 



alignments, is that star formation takes place in regions where 
the shock bounding the radio jet, has passed. 

Recent observations seem to support this idea. Deep spectra 
of the radio galaxy 4C41.17 at z = 3.8, show that the bright, 
spatially extended rest-frame UV continuum emission is unpo- 
larized and contains P Cygni-like absorption features, indicat- 
ing the presence of a large population of young, hot stars (Dey 
et al. |1997[ ). [Bicknell et al] ( 2000| ) argue that this can best be un- 
derstood if the shock associated with the radio jet has triggered 
star formation within the emission line clouds. 

A nearby example where stars might be formed under the 
influence of a radio source is the case of Cen A. Here, young 
stars are found near filaments of ionized gas in a radio lobe 



( Mould et al. 2000) 



Rees ( |1989| ) and [Begelman & Cioffi| ( |1989| ) analytically 
explored the evolution of intergalactic medium (IGM) clouds, 
overtaken by shocks from the cocoon of a radio jet. They argue 
that these clouds would be compressed and then gravitation- 



ally contract to form stars. However, [eke (1999) claimed that 



the destructive aspects of the interaction between the expand- 
ing cocoon and the clouds would dominate the evolution of the 
clouds. In his scenario the clouds evaporate and their material 
mixes into the jet cocoon. 
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Given the complexity of the interaction between the clouds 
and the jet cocoon, numerical studies are a good tool to in- 
vestigate this problem. Although the 'shock-cloud interaction' 
problem was studied numerically before, none of these studies 
addresses the effects of radiative cooling, important for inter- 
galactic clouds. Here we present new results of a numerical 
hydrodynamic study of the shock-cloud interaction problem, 
including the effects of radiative cooling. 

In Sect. 2 we describe the general problem of shock-cloud 
interaction and the application to IGM clouds. Section 3 deals 
with the numerical method, and Sect. 4 contains the results, 
which we further discuss in the fifth section. We sum up the 
conclusions in Sect. 6. 

2. Shock-Cloud Interactions 

Many numerical studies of single shock-cloud interactions 



have been carried out. Woodward (1976) being one of the first 



Various others followed, of which we will only mention two 



more recent studies: Klein et al. (1994), who provided a thor- 
ough analysis of the problem, and [Poludnenko et al. (2002), 
who studied the case of a shock running over a system of 
clouds; see these two papers for an overview of the literature. 
It is notable that in nearly all numerical studies to date, radia- 
tive cooling was either neglected or had little effect. For work 
con sidering the large sc ale e ffects of the pa s sage o f radio jets, 
see ISteffen et al.| ( |l997| ) and Reynolds et al.| ( |200l| ). 

The evolution of a single, non-cooling cloud, which is run 
over by a strong shock wave, consists of three phases. Initially, 
the shock runs over the cloud. The time scale for this is the 
shock passing time, tgp = 2i?ci/f shock, where Rd is the cloud 
radius, and z^shock the velocity of the passing shock. 
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The second phase is the compression phase, in which the 
cloud finds itself inside the high pressure cocoon. It is now un- 
derpressured compared to its environment, and shock waves 
start to travel into the cloud from all sides. This phase lasts for 
a time tec = ^ci/^^s,ci; the cloud crushing time, where Vg ^i is 
the velocity of the shock travelling into the cloud. For a strong 
shock this velocity is of order Us,ci = i^shock/-v/X' which x 
is the ratio of the rid to ricnv, the densities of the cloud and the 
environment, respectively; see [Klein et al.| ( |1994| ) for a better 
estimate. 

The third phase starts when the shocks travelling into the 
cloud, meet and interact. This produces a rarefaction wave trav- 
elling through the shocked cloud material. The cloud, which 
was compressed by the shock waves, now starts expanding 
again, and soon afterwards is destroyed and mixes in with the 
surrounding flow. This typically happens in a few cloud crush- 
ing times. 

2.1. Cloud properties 



Following Rees ( |1989| ), [Begelman & Cioffi] ( |1989D , and 
McCarthy ( 1993 ), we assume the undisturbed clouds to be the 
cooler and denser phase of an ionized two-phase IGM, of which 
the low density phase has a temperature of Tjg = 10^ K and a 
density of rijg = 10^^ cm^'^. Assuming pressure equilibrium 
between the two phases, a cloud temperature of 10'* K gives 
a density of rid = 10 cm^'^. We choose an initial radius of 
3 X 10^° cm 100 pc) and hence the cloud mass is 9.5 x 
10^ Mp,. Following the analysis of Cygnus A by Begelman & 



Cioffi (1989), we take the Mach number of the shock bound- 



ing the jet cocoon to be 10, yielding Wshock= 3500 km s^ 
(0.01c). With these parameters we obtain tgp ~ 5 x 10"* years, 
ice = 8 X 10^ years, and Wci,s = 120 km s^^. 

The cooling time can be estimated from tcoo\ = Cv^^^/pd 
(see e.g. Kahn 1976), where C is a constant depending on the 
cooling processes, with a value of 6.0 x 10~^^ g cm~^ s"* for 
a gas in coUisional ionization equilibrium at solar abundances. 
With the numbers above one finds tcooi = 2 x 10^ years. This 
is the shortest time scale thus far, showing that cooling will 
dominate the evolution of the shocked cloud. 

It is instructive to derive a condition for which cooling will 
dominate. Using the expressions for t^c and tcooh we find that 
the condition tec > lOicooi can be rewritten as 



Mci > X 

/ Wshock 

V lO^kms" 
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This shows that cooUng dominates for a large range of values 
for Mci, Wshock, X, and ric. For our values of Mc\ and ric, the 
shock velocity needs to be above 17 000 km s^^, or the density 
ratio X below 15, for cooling not to dominate the evolution. 

3. Numerical method 

The calculations were performed with a two-dimensional hy- 
drodynamics code based on the Roe solver method, an ap- 



1995). Second order accuracy was achieved with the superbee 
flux limiter, which was made less steep by lowering the co- 
efficients from 2.0 to 1.2; taking 1.0 would correspond to us- 
ing the minmod flux limiter, see Sect. 20.2 in l^aneyl ( |l998| ). 
Better two-dimensional behaviour was implemented by using 
the transverse waves method as described by LeVeque (1997). 

In order to include the effects o f cooling, we used a cooling 
curve ( Dalgarno & McCray 1972), which gives the cooling as 
function of temperature, for a low density plasma in collisional 
ionization equilibrium. This is a reasonable approximation of 
the real cooling processes of astrophysical gases. The radiative 
terms were implemented using operator splitting, where the ap- 
propriate radiative losses and gains were added as a separate 
source term every time step. The heating rate is proportional to 
the density, and was set so that for the initial conditions, heating 
and cooling in the cloud are balanced. 

In order to deal with short cooling times, we subdivided 
the time steps into smaller fractions of the order of the cooling 
time when applying the cooling. We imposed a minimum tem- 
perature of 10 K. This approximately corresponds to the cos- 
mic microwave background temperature at the redshifts we are 
considering. We did not follow the ionization state of the gas, 
but assumed the material to always be in collisional ionization 
equilibrium. 

The geometry of the grid was either cylindrical (i?, z), 
assuming cylindrical symmetry, or cartesian [x, y), assuming 
slab symmetry. The use of two different coordinate systems 
helps in understanding the true three-dimensional nature of the 
flow. Cylindrical coordinates are the proper choice as long as 
the flow pattern retains its large scale character, i.e. during the 
initial phase of the interaction. However, when the cloud starts 
to fragment, off-axis pieces are represented by ring-shaped 
structures. Furthermore, there is a strictly imposed symmetry 
axis at the centre of the cloud. In cartesian coordinates the ini- 
tial conditions do not describe a spheroid, but rather a cylin- 
der. On the other hand, the fragmentation is more properly fol- 
lowed, and no symmetry axis is imposed. 

We ran two simulations: in run A the shock wave interacted 
with a spherical cloud (with the parameters from Sect. 2.1) on 
cylindrical coordinates, and in run B with an elliptical cloud 
(with a semi-major axis of 100 pc, axis ratio 1.5, the major 
axis at an angle of 45° with respect to the incoming shock, and 
all other properties the same as in run A) on cartesian coordi- 
nates. Using an elliptical cloud, rather than a spherical cloud, 
further reduces the symmetry. For both runs the cell sizes were 
0.486 X 0.486 pc, using 800 x 1600 (A) and 1600 x 1600 (B) 
computational cells. 

4. Results of the simulations 

Figure |l] shows the logarithm of the density for run A at times 
0.79 X lO'' and 1.1 x 10^ years. Figure||shows the same for run 
B[^ The cocoon shock wave came from the left, and passed the 
entire cloud at i = 5 x 10^ years. At < = 0.79 x 10^ years 
the shock waves travelling into the cloud have just merged 



proximate Riemann solver (Roe 1981; Eulderink & Mellema 



Movies of the entire density evolution of the two runs are available 
with the electronic version of this letter. 
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Fig. 1. Colour plots of logj^Q of the number density (cm^'^) att = 0.79 x 10^ years (left), and t — 1.1 x 10^ years (right) for 
run A (cylindrical coordinates, spherical cloud). The white contour indicates the original cloud position. Only a fraction of the 
computational domain is shown. The coordinates are relative to the lower left corner of each frame. The shock wave came from 
the left. 
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Fig. 2. As Fig. 1 for run B (cartesian coordinates, elliptical cloud). 
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(compare with the estimate for tec in Sect. 2.1). In the non- 
cooling case this is followed by a re-expansion of the shocked 
cloud (due to the extra heating generated in the merging of the 
shocks), but here the excess energy is radiated away, and the 
merging of the front- and back-side shocks leads to the for- 
mation of a dense, cool, elongated, but fragmented structure 
('sheet') perpendicular to the flow direction in run A, and more 
parallel to the major axis orientation in run B. In both cases 
there is a concentration near the centre of the former cloud. 

The two righthand boxes of Figs. |l] and ^ show how this 
sheet fragments further. In run A, the imposed symmetries lead 
to an elongated concentration of material on the axis, which we 



measured to contain approximately 10% of the original cloud 
mass. The rest of the cloud material is compressed into dense 
structures, spread out over a volume which is 30% of the orig- 
inal cloud size (part of the outer contour of the original cloud 
boundary is indicated in Fig. |l]). 

In run B the cloud develops into an ensemble of dense small 
fragments, filling an area of approximately the same diameter 
as the original cloud. Without imposing symmetry, the largest 
and densest fragment is found near the centre of the ensemble. 
The integrated density of this largest fragment was measured 
to be some 30% of the integrated density of the original cloud. 
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In both runs, at the end of the simulation, less than a per- 
cent of the original cloud material has been mixed into the co- 
coon, showing that the evaporation process is slow, as is ex- 
pected for high density contrasts. The ensemble does spatially 
disperse since the velocities of the fragments range from 90 to 
500 km s^^, the leftmost fragments having the lowest veloci- 
ties. 



5. Discussion 

The simulations presented here show a comple tely n e w be - 
haviour compare d to the sce narios presented in Rees ( 1989 ) 
(compression), or Icke ( 1999 ) (disruption). In our simulations, 
instead of being simply compressed or disrupted, the cloud 
breaks up into many small dense fragments, spread out over a 
certain volume, and which evaporate only slowly. This has not 
been seen in numerical simulations before. It is completely due 
to the introduction of cooling, which prevents the 'third phase' 
or re-expansion (see Sect. 2). 

We expect that full three-dimensional simulations will 
show a result lying somewhat in between what we found in 
runs A and B. It would definitely enhance rather than suppress 
the fragmentation, since there is one extra degree of freedom 



available for instabilities (see e.g. Xu & Stone 1995) 



There are a number of processes which could work against 
the cooling, and hence slow down the compression. These 
are for example heating by the UV and X-ray photons from 
the AGN and the presence of a magnetic field in the clouds. 
Simulatio ns of magnetized flo ws in three dimensions, as re- 
ported by Gregori et al. (1999), show that if the magnetic field 
is strong enough, it will actually enhance the fragmentation of 
the cloud, and presumably aid evaporation rather than compres- 
sion. However, these simulations did not include the effects of 
cooling, so it is difficult to compare their results to ours. 

Note that whenever the cooling time is substantially shorter 
than the cloud crushing time, we expect an evolution similar 
to the one above. Equation 1 shows that this holds for a wide 
range of cloud paramete r s. For example, the interstellar clouds 
from [Poludnenko et al. (2002) should strongly cool, be com- 
pressed, and develop into a long-lived mass loading flow, some- 
thing which these authors failed to achieve in their non-cooling 
simulations, where the clouds are destroyed within a few tec- 

The further evolution of our fragments will be dominated 
by two processes: gravitational collapse, and further accelera- 
tion and erosion by the passing flow. All fragments found in 
our simulations will collapse under their own gravity, which 
makes them smaller, and even harder to disrupt and/or acceler- 
ate. As pointed out in Sect. 4, nearly all of the original cloud 
material ends up in these dense fragments, and would be avail- 
able for star formation. This implies that the estimate for the 
induced star formation rate from Begelman & Cioffi (1989b, 



is still valid. For a cloud filling factor (by volume) of 10 
and a relativistic jet, they find an induced star formation rate of 
~ 100 M0 yr^^, in rough agreement with the observations. 



6. Conciusions 

We have for the first time simulated the cooling dominated evo- 
lution of an intergalactic cloud which is overrun by the cocoon 
of a passing radio jet. Previous analytical studies conjectured 
that the cloud would either be compressed, or be completely 
destroyed and evaporate into the cocoon. We instead find a 
new picture. Radiative cooling is so rapid, that nearly all of the 
cloud mass is compressed into many small and dense fragments 
with a long hydrodynamical survival time. These fragments are 
likely to collapse and form stars, in line with the scenario of jet 
induced star formation. 

This type of fragmentation is expected whenever the 
cooling time is much shorter than the cloud crushing time. 
Evaluating this condition, shows this to be case for a wide 
range of parameters, stretching from intergalactic to interstel- 
lar conditions, see Eq. 1 . The collapse-and-fragment sequence 
we find, may well be the way to create long lived mass loa ding 
flows inside post-shock regions ( Hartquist & Dyson 1988). 

These simulations are only a first step, and definitely more 
work is needed. In future papers we plan to explore the effects 
three-dimensionality and self-gravity have on the fragmenta- 
tion process. 
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